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1. Introduction 

Many practically important control systems do not have a 
stable equilibrium. Instead, they may perform satisfactorily 
while possessing only the kind of behavior characteristic that 
has been described under various boundedness or practical 
stability definitions [1]. This report discusses a class of 
techniques which were developed for obtaining quantitative 
information about the boundedness properties of such systems. A 
system of particular interest in this work is the sampled-data 
control of satellite attitude with quantization. Such a system 
will be employed throughout as an example of the application of 
the techniques described. 

The report begins with an introduction of various relevant 
stability concepts as a series of definitions in Section 2 ; 
interrelationships are discussed between various definitions in 
common use. Also included here is a description of the model of 
a basic sampled-data control system with quantization. 

Section 3 describes the basic technique of estimating 
boundedness regions by means of quadratic Liapunov functions. It 
also states a sufficient condition, for a certain class of systems, 
for the existence of a boundedness region. The proof of this 
condition (Theorem 3-1) involves the estimation technique. Some 
of the results in this section (Eqn. (3-12) and Eqn. (3-14)) 
parallel results of Johnson and Lack [2,3], which have also been 
applied recently by Parker and Hess [4]. 
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Section 4 applies a particular quadratic Liapunov function 
to the Lur' e-Postnikov class of systems, along the lines proposed 
by Weissenberger and Siljak [5, 6], This class of systems is more 
restricted than that to which the techniques of Section 4 applies ; 
however, the stability property established is the stronger one 
of absolute boundedness (Definition 4). Relationships are 
developed for applying the technique to the specific system with 
quantization. 

Section 5 presents an example of the calculation of region-of- 
boundedness estimates. 

Section 6 discusses the linear analysis of the system used 
in the example and simulation results of the nonlinear system for 
use in comparison with results obtained in Section 5 by Liapunov 
techniques . 

Section 7 presents conclusions. 

The Appendix contains a paper*"On Practical Stability" , by 
Ljubomir Grujic, Visiting Research Associate in the Electrical 
Engineering Department, University of Santa Clara, on leave from 
the Mechanical Engineering Department, University of Belgrade, 
Belgrade, Yugoslavia. This work was motivated by an interest in 
applying the practical stability concept to the systems of this 
report, and although it treats only continuous -time systems, 
useful extensions to discrete-time systems are obvious. Connec- 
tions between boundedness and practical stability concepts are 

discussed in Section 2. 

35 

The paper was presented at the Fifth Asilomar Conference 
on Circuits and Systems, Pacific Grove, California, November 8-10, 

19 71. 
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Notation 

Throughout this report, except where otherwise noted, lower 
case Roman letters denote vectors , capital Roman letters denote 
matrices, lower case Greek letters denote scalars, and capital 
Greek letters denote sets (except the letter X, which represents 
the set of all points in the state space). Vectors will be 
considered as column matrices. The superscript T denotes the 
transpose and * denotes the conjugate transpose. The notation 
H > 0 means that H is positive definite real symmetric matrix. 

The letter t is used for discrete-time index, and the letter V 
for a Liapunov function. The region Q is the complement of ft. 
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2 . Preliminaries: Definitions and Models 

a. Definitions 

Consider a discrete-time system described by the equation 
x t+l = + f(t, x t ), t = t Q , t Q + 1, (2-1) 

where x^ is the n-vector state of the system and g and f are 
n-vector functions of time t and state x^ . The vector f is 
considered as the input to the fundamental, unforced system 

x^. + ^ = gCtjX^) C 2 — 2 ) 

In later applications in this work, a special form will be 
taken for the function g , 

gCt, x t ) = P x t + q K° r t >i cr t = r T x t » (2-3) 

a decomposition into a linear part and a special nonlinear part: 

P is an n x n matrix of constant coefficients, q is an n-vector of 
constant coefficients, and <J> is a scalar function of the linear 
combination cr of system states. However, for purposes of stating 
definitions in this section we for the most part retain the more 
general system description (2-1). 

Definition 1 

All the motions x -t-(x ,t ) of system (2-l) / are bounded , if 

for each initial state and time (x ,t ) there exists a number 

o o 

6(x .t ) > 0 such that 
o o 


|x t (x o’ t o )| * 6 * t - t o • 

The motions of system (2-1) are then said to be Lagrange s table . 

By itself, such a boundedness property may give little useful 
information about the behavior of the system; one quite often 
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desires at least that all motions ultimately satisfy a particular 
bound. We are then led to construct 
Definition 2 

The motions of system (2-1) are said to be ultimately 
bounded if there exists a number 6 such that for each (x , t ) 
there exists at, > t such that 


l x t Cx o ,t o ) I < 5 » t 1 t i • 

Thus, ultimate boundedness implies the existence of a 
bounded region ft containing the origin which all solutions 
ultimately enter. 

It is frequently useful to consider a modification of 
Definition 2 to explicitly recognize this region ft and also to 
take into account the fact that for some systems the ultimate 
boundedness property is not global with respect to initial states : 
that there are system states from which motions do not enter the 
region ft. 

Let ft jC ft ^ be bounded regions containing the origin. 


Definition 3 

The motions of system (2-1) are said to be ultimately bounded 
with respect to regions ft^ and ft 2 if for each t Q and for each x q e ft 2 
there is at, > t such that the motion x, Cx ,t ) e ft, for all t > t, . 

-L — O lO Ol — 1 

A modification of Definition 3 will find application later to 
the special system of the form (2-1) - (2-3). Suppose that $ is a 
certain class of nonlinear functions 4>(a^.) and that F is a certain 
class of input functions f(t, x^.). The class $ will be described 
later, the class F is the class of bounded inputs, 
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F: |f(t,x t )| < V . (2-4) 

Definition 4 

The motions of system (2-1) - (2-3) are said to be absolutely 
ultimately bounded with respect to regions and fig » and to "the 
class of nonlinearities <J> and to the class of inputs F, if for each 

t ,x, each x e fi 0 ,each 4>e$,and each feF, there is a t,> t such that 
o o 2 ’ lo 

the motion x. (x , t ) e Q, for all t > t, . 

tool — 1 

Closely related to our definitions of boundedness are several 
notions of practical stability , which were stated originally by 
LaSalle and Lefschetz, Let Q q CQ be closed and bounded regions 
containing the origin. 

Definition 5 

The motions of (2-1) are said to be practically stable if for 

each feF, each t , and each x^ e (D , the motion x. (x , t ) e Q for 

o o o to o 

all t >_ t Q . A stronger stability is described by 
Definition 6 

The motions of (2-1) are said to have strong practical stability 

if they are practically stable and if in addition for each feF, each 

t , and each x , there exists a t, > t such that the motion 
o’ o’ 1 o 

x t (x , t ) e Q for all t >_ t^. 

Figure 2-1 illustrates for comparison purposes three divisions of 
stability behavior: ultimate boundedness with respect to 5)^, > 

practical stability, and strong practical stability. 

Note that the roles of the region of initial conditions (J^* 

Q q ) and the region in which solutions ultimately enter and remain 
(fl^» Q) are reversed in ultimate boundedness and practical stability: 
in ultimate boundedness there is a kind of convergence toward a 




Strong Practical 
Note : The indicated behavic 

FIG. 2-] 
COMPARISON OF STAE 



Def. 1. Lagrange stabili^ 


Def. 2. Ultimate boundedness, with number 6 


Ultimate boundedness with respect to Q 2 , where = X, 

the entire state space " " 


Def. 3. Ultimate boundedness with respect to tip ft 2 , ^ 2 i X 


Def. 5. Practical stability with respect to Q , Q 


Def. 6. Strong practical stability with respect to Q q , Q 


Interrelationships Between Stability Types 
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neighborhood of the origin in that while in practical 

stability, motions are allowed to enter a larger region than the 
region of initial states, QOQ q . Strong practical stability, 
however, combines elements of both: there is both practical 

stability with respect to Q q , Q and ultimate boundedness with 
respect to Q , X, where X is the whole state space. 

Figure 2-2 shows interrelationships between various stability 
types. The arrows mean "implies" and are to be understood, for 
example, in the following sense: ultimate boundedness of motions 

with respect to certain regions Q 2 (ft 2 ^ X ) implies that there 

exists some other regions Q , Q for which the motions are 
practically stable, 
b. Models 

Systems of particular interest in this work are sampled-data 
control systems containing quantization nonlinearities. The block 
diagram of a class of such control systems is shown in Fig. 2-3; 
important features are: linear plant dynamics G(S) , time delay T^, 

sample and zero-order hold, and a single quantization nonlinearity, 
which is shown in detail in Fig. 2-4 in a form with saturation. 

In subsequent developments , it will be necessary to put 
system equations in the state -variable form (2-1). Define the 
discrete-time open-loop transfer function of the system of Fig. 2-3 
as 


where 


-2lc k } 

$(z) = ^<t> ok Ca(k)]} 


(2-5) 

( 2 - 6 ) 

(2-7) 


From Fig. 2-3 we calculate 


G(z) 



-Ts 


e 



G( s) } 


s 


( 2 - 8 ) 








1 ] 



FIGURE 2-4 

QUANTIZER NONLINEARITY WITH SATURATION 
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where 


Then 


= ( 1-z” 1 ) ^{e _T d S GC s ) } 
G(s) 


Gi(s) = s 


G(z) = (1-z ^ } l m= i_ Td / T 

where ^ is the modified z-transform. 

0m 

The computation of (2-11) will yield in general, 


G( z) 


n-1 
£ b . z 
i=o 1 


i-n 


n-1 

1 + £ a. z 

i=o 


l-n 


, a^ = a^(m), = b^(m) 


(2-9) 

( 2 - 10 ) 

( 2 - 11 ) 


( 2 - 12 ) 


From (2-5) and (2-12) we can construct the state equations 

x t+l = p x t + G <J> Ccr t ) 

T 

°t = r x t 


(2-13) 


where 


P = 


0 1 

. 0 


0 0 

■a -a, 
o 1 


. 0 

. 0 


. 1 

. -a 


n-1 


; q = 


o 

o 

0 

■ 1 


; r = 


n-1 


(2-14) 
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3. Quadratic Liapunov Function Estimates of Boundedness Regions 
In this section we establish a sufficient condition for the 
existence of a region of ultimate boundedness with respect to 
X. In the process of establishing this condition by means of 
quadratic Liapunov functions , we formulate techniques for calculating 
estimates of the region 

First, consider the basic 

Theorem 3-1 : Suppose that the system (2-1) can be put into the 

form 


x t+l = Px t + f(t ’ x t ) 


( 3-1) 


where P is a constant coefficient, Hurwitz matrix and the vector f 
is bounded, 


feF, V t > t , Vx. 

— o 

Then there exists a region such that the system (3-1) is 
ultimately bounded with respect to X. 

Proof : Consider the quadratic scalar function 

V(x) = x T Hx 

and its difference along motions of (3-1), 


AV = V(x t+1 ) - V(x t ) = (Px t + f) T H(Px t + f) 

AV = - x t T Qx t + f T Hf 
where Q = H - P T HP . 

Since by assumption 


(3-2) 

(3-3) 


|X i (P)| < 1, i = 1,2 ..., n, 

for each Q > 0 there exists a unique solution H > 0 to Eqn. (3-3) 
[ 7 ]. If Q > 0 is chosen, the first term in (3-2) is a positive 



definite quadratic form. Since feF, the second term is bounded. 
Hence , for sufficiently large |x| , AV < 0; in this region where 
AV < 0 all motions cross surfaces V(x) = constant from the outside 
to the inside. Consequently, an estimate ^ of the region will 
be given by 

= (x: V(x) < v 1 ) (3-4) 

n i rb ^ (3-5) 

where = maximum V(x) (3-6) 

x: AV(x,t) = 0, t > t . 

— o 

Note the conservativeness of the estimate (Eqn. (3-5)). 

This demonstration has established the theorem and also provided 
a procedure for calculating estimates 
Procedure : 

i) Choose Q > 0 

ii) Determine H from the equation 

H - P T HP = Q (3-6) 

iii) Calculate v^: 

= maximum V(x) 

x: AV(x,t) = 0, t > t ' (3-7) 

— o 

. . T 

i.e., v-^ = maximum x Hx subject to the constraint 
- x^Qx + f(t ,x)^Hf (t ,x) = 0, t _> t 

iv) = x: V(x) < v 1 (3-8) 

In all but contrived, low-order examples, this procedure 
requires a computer. The calculation of H in (3-5) is straight- 
forward: the equation may be transformed into a set of linear 

equations in the elements of H, or else various direct algorithms 
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for the solution of the continuous time version of (3-5) [8] may 
be employed in conjunction with a bilinear transformation [8]*. 

The calculation in step iii) contains the most potential 
difficulties, depending to some extent on the nature of the 
function f(t,x). References [9], [10], [11] discuss computer 
methods which were used successfully in solving the analogous 
equations which arise in computing quadratic estimates of regions 
of asymptotic stability for the time-invariant case, f = f(x) . 

The quality of the approximation, i.e., the "closeness" of 
to will depend in general on the choice of Q (as well, of 
course, on one's criterion for evaluating its quality); for each Q 
a different may be expected to result. A strategy for obtaining 
improved results would be to determine the elements of Q (subject 
to the constraint Q > 0) which extremize some measure of the quality 
of The volume is a reasonable quality measure and is readily 

calculated for a quadratic. The resulting modified procedure is 
analogous to that used by Weissenberger [12], Nelson [10] and 
Geiss, et al., [ll] to calculate estimates of regions of asymptotic 
stability.** Its defects are the likelihood of excessive computer 
time and convergence problems for high-order systems. 

Difficulties in the calculation of in (3-7) may be avoided 
by introducing some degree of approximation and accepting more 
conservative results. The first term in (3-2) satisfies the 
inequality 

- x^Qx < - Aq |x|2 (3-8) 

Let B = I + 2 (P-I) -1 and Y = (B T -I)Q(B-I) . Solve B T H + HB = - 

** 

In the case of boundedness regions, the volume would be minimized 
as opposed to the maximization in the case of regions of asymptotic 
stability. 


hOjKl 
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(3 


where X^ > 0 is the minimum eigenvalue of Q. Because feF , the 

second term satisfies the inequality 

T 2 

f Hf > X h y Z 

where X^ >0 is the maximum eigenvalue of H. 

Eqns . (3-2), (3-8) and (3-9) give 

AV < - X Q |x| 2 + X h y 2 


(3 


AV < 0 


for all x such that 


or 


- Xq } x | 2 + X h y 2 < 0 

I x | >Y^ /2 Xq 1/2 = p 


Consequently, an estimate C S C fh is given by 


( 3 


where 


8^ = (x: V(x) < v x > 


v, = maximum V(x), 


|x|= p 

The number is readily found to be given by 

• , 2 

V 1 = V 


or 


ss 2 2-1 

v i = y 


(3 

(3 


(3 


In situations where the set of initial states is not 
whole state space X, we make use of 
Theorem 3-2 ; 

Let V(x) > 0, and define the (bounded) regions and 


and 


^ = (x: V(x) < v } 

&2 = ^ x: V(x) < v 2 ) , 
0 < v -^ < v 2 . Ass ume 

^ 2 1 ^ x: < 


-9) 

- 10 ) 

- 11 ) 

- 12 ) 

-13) 

-14) 

the 

as 


where 
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Then, the motions of the system (3-1) are ultimately bounded with 
respect to the regions and ^ • 

Note that this procedure will produce results which are 
conservative in the sense that 

and ^2 CT ^2 , 

Theorem 3-2 will be employed in Section 4. 
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4 . Lur * e-Postnikov Liapunov Function Estimate of Boundedness 
Regions 

In certain applications , attitude control systems may be 
modeled by the Lur ' e-Postnikov class of systems where the linear 
part of the system is not asymptotically stable and the quantizer 
represents the nonlinear characteristic. A simple transformation 
can be used to make the linear part asymptotically stable which, 
in turn, forces the quantizer characteristic to violate the usual 
sector condition in the neighborhood of the origin. An approach 
to the analysis of this class of systems is to estimate the result- 
ing regions of ultimate boundedness as proposed in references [5] 
and [6]. 

The estimation procedure makes use of a quadratic Liapunov 
function, a modification of the Tsypkin frequency criterion [1.3], 
or algebraic test Il4],and the Szego-Kalman construction [15]. In 
applying the procedure to a specific situation, one has several 
parameters available with which to improve the estimates of the 
regions of boundedness. 

We consider a free, discrete-time system of the Lur'e class 
described by the nth-order difference equation, 

x t+l = Vt + q W + 

a t = r x., t = 0,1, ... ( 4-1) 

where x, , q, and r are real n-vectors ; P is a real n x n matrix; 
t o 

4> 0 (a) is a real scalar function of the real variable a, which may 
have isolated discontinuities; and f(t,x t ) eF. It is assumed that 
the pair (P Q ,q) is completely controllable and the pair (P Q ,r ) is 
completely observable. 
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Note that in this section, unlike the previous one, we identify 
a particular scalar nonlinearity (although we retain also a vector 
nonlinearity f ) . 

The system (4-1) is transformed into 


x t+l = Px t + + 

T 

a = r x^. , t = 0 ,1 , . . . 


(4-2) 


where 


T 

P = P + Tqr 
o 


(4-3) 


<f>(a) = <j> (a) - xa 
o 

The number t in (4-3) is chosen in order to insure that the trans- 
formed matrix P is Hurwitz, that is, the n eigenvalues A k (P) of 
the matrix P satisfy 


| X R (P) | < 1, k = 1,2 , . . . , n (4-4) 

and to guarantee that the transformed nonlinear function <(>(a) 
belongs to the class $ defined by 

$: 0 < o<t>( o) < ko , £|a| < a 2 

| <fi(a) | ^ g , |a|<a ]L ( 4-5 ) 

where a^, g >. 0 , 0 < 02 £ + ” , and the numbers k > 0, 6 > 0 are 
selected to satisfy the inequality 


-1 


+ Re “6 h*h > 0, Vz:|z 


(4-6) 


In (4-6), 


X(z) = r T (P - zl)" 1 q (4-7) 

is the open-loop transfer function of the linear part of (4-2) 
from the input 4>(a^.) to the output - a , and h(z) is the complex 
vector defined by 

h(z) = (P - zl ) -1 q (4-8) 

Condition (4-6) is necessary and sufficient [16] for the existence 
of a function 
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V(x) = x T Hx (4-9) 

with H > 0, such that along the solutions of (4-7) 

- AV = 6 1 x | ^ + (£<J> + P T x ) 2 + (p(o - k ^ <j>) 

- 2x T Hf (4-10) 

where the matrix H, vector P, and scalar 5 satisfy the equations 

H - P T HP = PP T + 61 

q T Hq + ? 2 = ic ” 1 (4-11) 

P T Hq + qc - j r 

Since the system is described within a Lur'e context, it is 
natural to use a boundedness property which reflects the class $ 
as well as the class F. Such a property was characterized in 
Section 2 in Definition 4 as absolute boundedness with respect to 
regions ^ and an d classes $ and F. To calculate estimates 


and ^2 °f the regions and we ma ke use of the function 
V(x) and Theorem 3-2. Use is also made of the following 
inequalities, obtained from (4-10): 


- AV > 6 | x | 

- A V > 6 I x I 


- I x I - y 

- * h Y I x| , 


a | < a 
a | > a 


(4-12) 


where X^ is the largest eigenvalue of H, and 


y = 0(a 1 + k " 1 &) (4-13) 

Based upon inequalities (4-12), the least conservative estimates 
obtainable for system (4-2) with properties (4-4) and (4-5) are 
the following: 


A A 


(2^ - 

(x: 

V(x) 


A 

^ 2 = 

fx: 

V(x) 

A 

< v 2 > 


(4-14) 
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= max x Hx 
xe A ^ 

A fji 

v 0 = min x 1 Hx , 


(4-15) 


A 1 = r i U r 2 U r 3 

r i = (x: I X I = P 1 , |r T x| < c^} 

r 2 = {x: I X I = P 2 , | r T x | > a x > 

r 3 = {x: | r x | = a-p p 2 < |x| < p x > 

X H y . r , X H Y .2 j. 4-.1/2 

P 1 " TT + S ] 


(4-16) 


(4-17) 


P 2 = T" 


A 2 = {x: |r x| = a 2 > 


(4-18) 


An illustration of the sets A 1 and A 2 in two dimensions is given 
in Fig. 4-1. 

A 

The value of v^ given in (4-15) can be calculated by the following 
procedure: 


(i) Define 


where 


v, . = max x Hx , 


A, . = {x: | x | = p, } 


(4-19) 


(4-20) 


From (4-19) and (4-20), we compute 


V 11 = X H P 1 


Using (4-21) we test 


"1/2 T >.-1/2 T 
( A^e e ' r e 


> o , v 1 -- V U 


< 0 , go to (ii) 


where e is an eigenvector corresponding to A. 


(4-21) 
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(ii) Proceeding from (4-2 2b) we calculate as 
follows : 


12 


p l e i + 2 ct x e 2 


where 0^ and 0 2 are solutions of the equations 


Hx - 0^6x - j 0 2 r = 0 


T 2.-1 

X X = p^o 

T 

r x = a. 


(4-23) 


Also 


Then 


V 13 = X H P 2 * 


v 


11 


= max {v 12 , v n q } . 


13‘ 


(4-24) 

(4-25) 


Note that the value can always serve as a conservative 

' A 

approximation for v^. 

The value of v 2 given in (4-15) is obtained simply as 

(4-26) 

Note that a necessary condition for the existence of estimates 

A A . . 

^ 2 , and a sufficient condition for the existence of regions 


v 2 = a 2 (r T H" 1 r)“ 1 . 


^ 2 is that 


V 1 < v 2 


(4-27) 


It is of interest to point out certain special cases of this 
analysis that arise when the constraints on nonlinearity and a 2 
assume certain limiting values. If is reduced to zero the above 
procedure produces finite regions of absolute stability. When a 2 
is infinite the procedure establishes either a global property of 
absolute ultimate boundedness (ci^ 0), or a global property of 
absolute stability (a^= 0). 
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Application to the Quantization Nonlinearity 

In order to apply the foregoing analysis technique to the 
computer control system with quantization described in an earlier 
section (see Fig. 2.3 and Fig. 2.4), we require calculation of 
the numbers a^, , and 3* Figure 4-2 shows how these quantities 

depend on the transformation parameter t and the Popov sector 
parameter k. 

Inspection of Fig. 4-2 shows that a-^ is determined by 
either a lower or upper sector intersection with the transformed 
nonlinearity , 

= max {a L , a y } , (4-28) 

where a L arises from the intersection with the lower limit of the 
sector (the a axis) and arises from the intersection with the 
upper limit of the sector (the line kct). The quantity a L is 
obtained from the relations 

a L e (4-29) 

for that integer value of n such that 


2n 


2n+l - T - 2n+3 
The quantity a„ is obtained from the relations 


2n+2 


a T 


n 


(4-30) 


(4-31) 


U K+T 

for that integer value of n such that 

2n-l „ n 2n+l 
~2 - ic+r - ~~2 ' 

A simpler, conservative estimate of may be obtained 

using the straight line envelope of the nonlinearity (the parallel 


(4-32) 



FIGURE 4-2 


TRANSFORMED QUANTIZATION NONLINEARITY 
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dashed lines in Fig. 4-2). We still use (4-28) but in place of 
( 4—29 )— ( 4—32) we have 


a L = = 271 '-t) 

and “U = a U = 2 ( T+T-l) 

Since it is necessary that > 0 , we have from (4-33) and 
the following limitation on the transformation parameter x: 


(4-33) 

(4-34) 

(4-34) 


1-T < X < 1 


(4-35) 


The approximate analysis from the straight-line envelope 
also gives, for 6, 

B = (l-x) a x + | (4-36) 

The value of a 2 is given simply by 

r 


a. 


m + 


x 

+ oo 


, x > 0 

, x < 0 


(4-37) 


In those applications, such as satellite attitude control, 
the linear part of the system is not asymptotically stable, the 
calculation of regions of absolute ultimate boundedness requires a 
transformation described in (4-3). Consequently, the first step in 
the calculation is a selection of x such that P is Hurwitz and (4-35) 
is satisfied. To verify (4-35)we need a k which satisfies 
frequency condition (4-6) for a certain choice of 6. Once k and 
6 are selected, one calculates the vector g by the Szego-Kalman 
[15] factorization procedure and computes the matrix H from Eqs . 
(4-11). After H is determined, one calculates the numbers v 1 and 
\>2 as explained at the end of the preceding section. If either 
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V 1 > v 2 ’ or es ‘ t ^ mates and ^2 are nat satisfactory, one repeats 
the entire procedure with a different choice of the transformation 
t, and possibly different numbers for 6 and ic. As is clear from 
this outline, the application of the proposed method to higher 
order systems would require utilization of a computer. 
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5 . Example 

In this section we perform a numerical calculation of a 
boundedness region by Liapunov methods. The result will be 
compared with that to be obtained in Section 6 by simulation. 

Let the specific system be that to be described subsequently in 
Section 6 in Fig. 6-7a; this model can be taken as a very 
simplified representation of a satellite attitude control system. 
Direct calculation (by integration between sampling intervals) 
gives the state equations 


where 


x 


t*l = P o x t + =1 W 


T 

r x. 



T] 

GK 1 ] 


(5-1) 


and <J> is the quantization nonlinearity of Fig. 2-4 with m = °°. 

The quantity BETA in Fig. 6-7 has been given a unity value, and 
the quantity GAMMA is denoted above by G. 

Using the methods of Section 3 we require a transformation 
of Eqn. (5-1) to form in which the nonlinearity is bounded. Such 
a transformation is given by Eqn. (4-3) with t = 1. Equation (5-1) 
then becomes 

x t + l = P x t + q cf>Ca t ) 


(5-2) 
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where 


P = P + q r 
o H 


and 


<p = d> - r x 
o 


To calculate a boundedness estimate let us use the simplified 
form of the calculation which relies on Eqns . (3-6), (3-12), and 
( 3-14) : 


H - P HP = Q > 0 
H = ^x: V(x) <v^* 


(5-3) 

(5-4) 


2 , 2.-1 

Y X h Xq 


(5-5) 

where V(x) = x^Hx (5-6) 

The specific choice of numbers in the simulation in Section 

6 , 


gives 


K = 

0.292 


0 

It 

i — 1 
* 

1.146 


G = 

1.6 


T = 

1.0 






0.76 

0.08 

P = 




• 

o 

1 

CO 

CO 

• 

o 

1 


L 


q T = [-0.5 

-1.0] 


r = [0.467 
Y = 0 . 56e 


1.83] 


Choosing 


Q = i 


and solving (5-3) , we obtain 


H = 


2.52 

0.78 


0.78 

2.87 
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Eqn. (5-5) then gives 

« 2 
= 3 . 8e 

so that the estimate of the boundedness region is given by 

x < 3. 8e^ (5-8) 

This region is shown in Fig. 5-1 together with the region as 
estimated from the simulation results of Section 6 in Figs. (6-13) 
through (6-15). The overlap of the regions may be due to inaccura- 
cies in estimation from the simulation results : the actual region 

may be smaller than was concluded from the simulation, due to the. 
very complex, long-duration dynamic behavior of the simulated 
system. In any case, more examples and comparisons are needed 
before a definitive statement can be made about the merits of the 
Liapunov approach. 


2.52 
0. 78 


0. 78 
2. 87 
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6 . Linear Analysis and Simulation Results 
a. Linear Sampled Data Studies 

One model of an attitude control for a satellite has been 
given by Seltzer [17] and is represented in Fig. 6-1 as Model A. 

An alternate representation representing somewhat different 
instrumentation is given in Fig. 6-1 as Model B. Both models are 
linear sampled data systems with digital compensator and a time 
delay due to the use of the digital computer. It is intended that 
a nonlinearity (a quantizer) be inserted into Model B, and the 
effect of such nonlinearity on stability is to be studied. First, 
however, a linear analysis of Model B is undertaken with appropriate 
comparisons to Model A. 

Symbol equivalences in the two models: 

Model A Model B 


K, 


K 


K. 


K 


K, 


K 


K 


k Q = K q TV2I 


k Q = KKqT / 21 


k 1 = K 1 T/I 


k ± = KK y T/I 


In Z-domain, the characteristic equations and transfer func- 


tions for Model A (with digital filter and F(Z) = 
B (without digital computer) are as follows: 

Model A 

Characteristic Equation: 


T+K, 


-) and Model 


Z 3 + (K 2 - 2 ) Z 2 + (1 + k Q + k x - 2K 2 )Z + (k Q - k ± + K ? ) = 0 
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MODEL A (Seltzer's model) 


Vehicle 

-X- 

Digital 

Control 

Law 

T 

Filter 

V K 1 S 


F(Z) 




Zero 

Order 

Hold 

- 

Computa- 

tional 

Delay 

- 

, -si 

1 - e 




S 


e- T D S 



Vehicle 

Dynamics 


IS 


MODEL B 



FIG. 6.1: TWO MODELS FOR A LINEAR SYSTEM 
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Transfer Function: 

<k Q + k x )Z + (k Q - k 1 ) 

Z 3 + (K 2 - 2)Z 2 + (1 + k 0 + k x - 2K 2 )Z + (kg - k 1 + Kg) 

(k Q + k 1 )Z" 2 + (kg - k 1 )Z“ 3 

1 + (K 2 - 2 )Z -1 + (1 + kg + k 1 - 2K 2 )Z" 2 + (kg - k x + Kg)Z 
Model B 

Characteristic Equation: 

Z 3 + (K 2 - 2 ) Z 2 + (1 + kg + k 1 - 2K 2 )Z + (k Q - + kg) =0 

Transfer Function: 

* _ ( 0 . 5 )KT 2 (Z + 1) 

Z 3 + K 2 - 2 ) Z 2 + (1 + kg + k x - 2K 2 )Z + (kg - k x + K 2 ) 

or (set K = 1) 

_ (0.5)T 2 (Z“ 2 + Z~ 3 ) 

1 + (K 2 - 2 ) Z -1 + (1 + kg + k 1 - 2Kg)Z -2 + (kg - k + Kg)Z 

Notice that both Model A and Model B have identical characteristic 
equations . 

In Z-domain, the stability boundary is the unit circle, i.e., 

| Z | = R = 1. 


R 


or : 


1 _ = 


Stability : | Z | < 1 

Instability: |z| > 1 


Z-plane 


Z = e ST = RE 10 


t, -CW T 
R = e n 


0 = W t(i-?- 2 ) 1/2 


5 
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where C = damping ratio 

W n = natural frequency 

T = Sampling period 

The stability contours (for complex boundaries) of K 2 (third 
parameter) in parameter kg - k^ plane are as shown in Fig. 6-2. 

This corresponds to Seltzer's Fig. 4-1, but is explicit, i.e., it 
has been computed quantitatively for a number of values, of K 2 * 

The plot is done by a digital computer by introducing Chebyshev 
functions into the characteristic equation. One real root 
boundary, Z = 1, is kg = 1 (i.e., K 2 - k^ plane); it is independent 

of the values of K 2 and in 3-D space. Another (Z = - 1) is 

calculated by setting Z = - 1 in the characteristic equation, and 
it is a plane defined by the equation k^ = 2(K 2 - 1) in the 3-D 
space which we can see is independent of k Q . 3-D space is not easy 
to show. However, in the parameter plane (2-D) for K 2 - = 1, the 
stable region is bounded by the two axes and the complex boundary 
as depicted in Fig. 6-2. Relative stability contours in the 
parameter plane K 2 = 0, K 2 = 0.5, K 2 = 1, K 2 =2 are computer 
plotted, shown in Figs. 6- 3, 4, 5 and 6, respectively, for discrete 

varying values of R and 6, again Chebyshev functions are introduced 

into the characteristic equation for the computer program. 

Figures 6-3, 5, and 6 correspond to Seltzer's Figs. 5,6 and 7, 
respectively, the only difference being in the variables chosen 
for mapping. 

The correlations of the above-mentioned graphs and those of 
Figs, -5, 6, and 7 of Seltzer's paper are developed as follows: 














Si/. 
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FIG. 6-3: RELATIVE STABILITY CONTOURS 
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O* T) " CW T , q t 7 m / n y. 2 \ 1 1 '2 

Since R = e n and 0 = W T(1 - C ) 

n 

log e R 2 1/2 

i.e., log e R = - ?W n T and — -- (1 - r) x/ = - ? 

i°§ p R o _i /o ? 

set — g — = A , we get £ = A(1 + A z ) W n T = 0(1+A ) 

Figure 6-7 is the computer generated data conversion chart to trans 

late between the R, 0, and £, W^T variables for db sign purpose. 

For example, choose K 2 - 1, R = 0.93, 0 = 1.2 in Fig. D, we 

find kg = 0.4, = 1.7, and in Fig. F we find ? = 0.07 and 

W^T =1.2. To check this in Fig. 6 of Seltzer's paper for kg = 0.4 

k^ = 1.7, one finds £ = 0.07 and W^T = 1.2. To choose a point 

within the unit circle of the Z-plane and measure R and 0 

(0 _< 0 < 2tt) we can immediately define its corresponding £ and 

W^T values for a certain third parameter K 2 . The process can be 

reversed, i.e., for a given £ and W^T , we can, through using 

Fig. 6-7, find the corresponding R and 0. The purpose of using 

R - 0 variables instead of ? - W T is for the convenience of 

n 

choosing a point in the Z-plane. Once W n T is found, for a 
predicted system natural frequency W , one can find the sampling 


1/2 


time T. 



ZETA 



FIG. 6.7: DATA CONVERSION CHART BETWEEN (R-0) AND U-W M T) CURVES 



b. 


Simulation Results 


One purpose of this project was to explore the possibility 
of applying stability theory to the definition of regions in the 
state space which the terminal state is guaranteed to reach from 
some defined set of initial conditions. Digital simulation was 
undertaken to provide "experimental" data which could be used to 
verify the conclusions drawn from the theoretical studies. This 
section describes the simulation program and some of the pertinent 
results . 

The block diagram used to represent the system in the 
simulation studies is given in Fig. 6.7a. This was modeled in the 
IBM 360 computer using the CSMP-360 program. The original program 
is given in Fig. 6-8, with a set of parameter numbers and initial 
conditions. Note that the two blocks, BETA and GAMMA, shown in 

Fig. 6- 7 a , were inserted for gain adjustment and distribution studies. 
It was thought that the location of gain with respect to the non- 
linear element might alter performance, but this was found to be 
untrue as far as stability is concerned. 

The numerical values selected, K0 = .292, K1 = 1.146, T = 1. 
and DEL3 = .1 were used rather than those supplied for the Sky lab 
because simple numbers were desired for the initial theoretical 
studies (which would include some long-hand calculations). Note 
that trapezoidal integration was used to avoid problems at the 
discontinuities, and the integration interval was DELT = .01, which 
turned out to be too large an interval. Fig. 6-9 a, b, c, show 
phase plane plots of system response to initial conditions. The 
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****CCNTINUCL'S. SYSTEM MQCCL ING PROGRAM**** 
***pppR|_E w INPUT STATEMENTS*** 


LABEL STAB! LI TV C MG 

PARAMETER KO=.292 - - — • 

PARAMETER Kl=1.146 
PARAMETER Bh TA = 1 . 

parameter gamma=i. 

PARAMETER T=l. 

-_P_ARAM.£ T E-R -GEL 3 = . 1 

PARAMETER P1=0. 

INCON XA = 1 • » >P = 2. 

HISTRY ZWOLC(ICC) 

xT=TNTGRL ( X A , X2 ) 

X2= INTGRL ( XB,BF;TA*A ) 

A~2 HCL C ( S» B-) 

B = CMZR(0LL3»C) 

C=GAMMA*CC 

C-C = ALPEA*D 

S = A L PHA 

ALP HA = I M R L L S ( P 1, TJ 

Z = M 1* X 2 

*=ir.c*xi 

METHCD TPAPZ 
PPSPAR XI. X2 

P-P-HW — X-l-»- X 2 » A — 

TIMER FIMTIM = 50i). ,CELT=.01,PRDEL = 1 . 
END 

STOP — 


FIG. 6-8: INITIAL PROGRAM 




FIG. 6-9 c: BETA = 1.0, GAMMA 








terminal portion of the trajectories shows what appears to be a 
limit cycle. The plots, however, are rather crude because the time 
increment was too large and the scale plots poorly chosen. It was 
decided to get plots with better portrayal of the cycles, so the 
program was rerun with better plot scales and with smaller time 
increment and smaller integration interval . Typical results are 
shown in Fig. 6-10. Note that the supposed limit cycles have 
disappeared. The obvious cause of the discrepancy is the integra- 
tion interval. This is mentioned in the report because the proposed 
Sky lab will use an on-line time shared digital computer, and the 
position and velocity measurements are quantized for use in this 
computer. Thus, if the integration interval or measurement 
granularity are chosen too large, a situation such as we have 
observed may be encountered. 

To continue the simulation studies the loop gain was varied 
to find a value for which the system would exhibit oscillatory 
characteristics. This was done to inverse the probability of 
finding limit cycles , since a part of the theoretical study is 
concerned with such phenomena. It was found that for relatively 
modest gain increases the system damping changes substantially. 

For K0 = .292, K1 = 1.146, BETA = GAMMA = 1, the system is 
heavily damped, exhibiting characteristics similar to the chatter 
mod- of relay servos. Changing only GAMMA, it was found that at 
GAMMA =1.8 the system appears to be divergingly unstable. It 
was decided that GAMMA = 1.6 was a suitable value for our 


purposes . 



FIG. 6-10a : 


BETA = 1.0, GAMI 



0.2 


















so 
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All remaining simulation runs were made with the forward 
gain set at GAMMA = 1.6, and the initial conditions were chosen 
on a square grid surrounding the origin of the phase plane. A 
copy of the program is given in Fig. 6.11; it differs from that in 
Fig. 6-8 only in the value of GAMMA, the integration interval, 
and the run time, print interval, etc. The results for this 
increased gain condition not only differ amazingly from those in 
Fig. 6-10, but showed some surprising symmetries which can be 
described briefly. First, the responses showed symmetry in a 
polar sense, i.e., initial conditions in the first and third 
quadrants gave responses which were identical in a polar symmetry 
sense, as were responses to initial conditions in the second and 
fourth quadrants. Initial conditions of position only gave 
comparatively well damped responses, with no chatter or limit 
cycles for small values of initial condition, but exhibited a 
chatter mode type of limit cycle* as it approached but did not 
reach the origin. Initial conditions of velocity only went 
immediately into a limit cycle type of oscillation about the 
origin, and as the magnitude of this initial condition increased 
several such modes of different amplitudes appeared. Initial 
conditions in the first and third quadrants gave limit cycles that 
did not enclose the origin, while initial conditions in the second 
and fourth quadrants gave limit cycles that did enclose the origin 

Because of the symmetries noted above, the results presented 
do not contain all of the data obtained but just representative 

“Note tne term limit cycle is used rather loosely here to 
describe a type of motion which is not precisely a limit cycle. 
This is discussed later. 



.**:?* CON T I NUCUS SYSTEM MODELING PROGRAM*** 
***FR0PLPM INPUT STATEMENTS*** 


LABEL STABILITY CMG 

PARAMETER KO* .292 

PARAVCT^P X 1 = 1 , 1 46 
PARAMETER B 0 TA = 1 . ,G AM YA = 1 . 6 
- P A R A Y rTFP T = I V •' 

PARAMETER D2L3=.l 
PARAMETER EO-s.il* 

INCON XA=0. »XP=-1 .6 
HISTRY — Z-H0L-B 


X 1 = I N T G R L ( X A , X 2 ) 

X2=TfTT GRL { XB rB ETft» A T~ 


A= ZHCL C ( ' 

>, RJ 

- BaCMZRil 

3£L3»£-> 


c=gayma*cc 

CC= AL PHA*C 

S**tPHA .. ... 

ALPHA =T YFUL S ( PI » T ) 

0= - Y - Z 

r=irr*xz 

Y=KC*X1 

■MET HC£l TRAP Z 

PREPAR XI, X2 
PRINT X 1 , X2 > A 

T^IMER — fW I M «6 9 > , D E L T « ,001, P R B EL ».l 
END 

STOP 


FIG. 6-11: REVISED PROGRAM FOR RESULTS OF FIGURES 6-12 

THROUGH 6-15 
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samples. Fig. 6-12 shows four trajectories starting from an 
initial position with zero initial velocity. Note that for small 
values of initial position there is no tendency to cycle, but at 
large values of initial position some cycling is produced, but 
does not enclose the origin. Fig. 6-13 shows three trajectories 
starting from an initial velocity with zero initial position. For 
a small magnitude of initial conditions, the system immediately 
starts cycling about the origin. For a substantially larger 
initial velocity two modes of cycling occur as the trajectory 
approaches the origin, but for a still larger initial velocity 
one of these modes disappears. Fig. 6-14 shows four trajectories 
starting in the first and third quadrants. All seem to terminate 
in a type of limit cycle which does not enclose the origin. 

Fig. 6-15 shows four trajectories starting in the second and 
fourth quadrants. In this case all seem to terminate in a type 
of limit cycle, but for the smaller initial conditions the cycle 
does not enclose the origin, while for the larger initial conditions 
it does. 

c. Comments on Simulation Results 

From the trajectories shown on Fig. 6-12 through Fig. 6-15, it 
is clear that the oscillations observed are not limit cycles in 
the usual sense, since the trajectory is not repeated exactly on 
successive cycles. There appears to be a "drift", i.e., successive 
cycles tend to be displaced along the position axis , usually 
tending toward the origin. There was evidence, in some of the 
print-out data, that the trajectories eventually reached the 
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FIG. 6-12b: BETA = 1.0, GAMMA = 1.6, XA = 0.4, XB = 0.0 
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FIG. 6-14a: BETA = 1.0, GAMMA = 1.6, XA= - 0.2, XB = - 0.2 
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FIG. 6-15d: BETA = 1.0, GAMMA = 1.6, XA = 3.0, XB = - 3.0 
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dead zone of the quantizer and that a no-power, constant velocity 
mode followed the cycling mode. Since the model has no damping 
in the dead zone, it was of interest to investigate whether the 
motion entered a very slow limit cycle mode crossing the entire 
dead zone or whether the trajectory entered a type of limit cycle 
around one edge of the dead zone. 

To explore these possibilities several of the preceding 
simulations were repeated using a longer problem time. These 
results were not conclusive, but several of the test cases did 
recycle about one edge of the dead zone without ever crossing 
the dead zone. 

No valid conclusions can be drawn at this point, especially 
in view of the fact that the simulation model was a much simpler 
system than any practical realization. It is clear, however, that 
the type of behavior observed is due to the sampled nature of the 
nonlinear system. One suspects that this oscillatory behavior 
can be changed substantially by altering the sampling period. 
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7 . Conclusions 

This report has considered the problem of estimating regions 
of boundedness for discrete-time dynamic systems. Based on 
Liapunov- funct ions , several methods were developed for this 
purpose. A technique based on simple quadratic Liapunov functions 
led to a number of possible variants , with various degrees of 
complexity and a wide range of numerical difficulty. An example 
by this method was performed and the estimate compared with one 
obtained by simulation. Tentative conclusions from this example 
are that Liapunov results may be good and that simulation results 
may be difficult to interpret and time-consuming to generate; more 
examples, however, will be required for a definitive judgment on 
the effectiveness of these methods. 

The other Liapunov-based technique made use of the Lur'e- 
Postnikov quadratic Liapunov function and yields estimates of 
regions of absolute boundedness. These results contain new and 
useful information regarding the influence of the nonlinearity on 
the boundedness region; this new information is apparently obtained 
at the cost of greater analytical complexity. The implications of 
this complexity, however, cannot be judged until further experience 
is gained with a computer implementation of this technique. 

It should be noted that this Lur'e-based technique is 
limited in its present form to a system containing a single non- 
linearity, while the simple quadratic Liapunov technique is 
applicable to systems containing multiple nonlinearities. 

Based on studies to this time, the Liapunov methods of this 
report remain potentially attractive compared to the simulation 
method. Further work would profitably be devoted toward 



71 


increasing experience with the application of these techniques 
to a specific system, and to the development of effective 
computer programs for this purpose. 
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APPENDIX 


ON PRACTICE, STABILITY 

• 1 

LJubosnlr T. :Gruj1d 
Mechanical Englhqpijlng Department 
University! of [Belgrade 
Belgrade,! aigoslavla 


ABSTWMjT . 

In this paper, a c. la** o& nonlinean. nonautonomou* system* with mu Itipte non- i 
linealities is c onsidened. SuHicient condition* one developed ion. a type of, 1 
pnactical stability with ipeciited iettlin# turn. The. condition* one. independent '■ 
o( the actual ionm of nonlinean chanactenisties so that the# can be intenpneted a* 
condition s ion. "absolute" practical & lability. The liability teat it reduced to j 
veniiication oi the Hunwitz pnopenty oi i a domtant matnix. T hi* make * the *tabil-. 
it y analysis ojj the c ontidened cla** ojj: iwrUinean systems convenient ion. machine ! 
computation*. ; f ] 1 

The pnopo*ed pnactical * labili ty analysis is applied to a tiUnd-onden *y*tem 
with ievenal noniineanitie* . 1 3 ! 


I 


1. INTRODUCTION 


((finite or Infinite) tine Interval. 


This paper Is concerned with a practical stability i 
analysis of nonautonomous nonlinear dynamic systems! \ ! 
which have a nunfcer of nonllnearltles depending on ; “ A 
several variables. 


2. NOTATIONS 

(a 1 j) , n « n constant matrix. 


Important results In practical stability analysis • 
over a finite time Interval have been obtained by 
Weiss and Infante [1] and Weiss [2]. They derived , 
necessary [2] and sufficient [1] conditions ex- 
pressed In terms of real valued functions. In this! 
work, as In reference [3], the function V(x)« b T |x|| 

Is proposed as a candidate for system Liapunov func- \ 'jH C R n * , H ■ {A: Re X.(A) < 0 , 1 ■ l, 

w j - * 

tlon. This function leads to algebraic conditions ^ ^, R n x „ j + ph* H • (h ) 
for testing pnactical sta bil ity with a ipediied ' * 


| jb • (b 1 b 2 ... b n ) T , constant positive vector. 

| ic * (cj c 2 ... c„) T , constant positive! vector. 

] V € R®, V ■ (z:| |z| | <_y) , set of all 
!j : tor disturbances. 

| |f:R n * R® ■ T ♦ R n , f ■ (f, f 2 ... f n ) 

I jg:* 1 * R® * T - R n , g - (g 1 g 2 ... g n ) r 


[allowed vec- 

(r 


t. 


, n). 


settling time., which Is defined over a prescribed 


I f N:R n ■* R n , N 


(N, N 2 


... N„) t . 


P(.)<=R n , ■ (x:|xi « s (,) ) , (•)! ■ a, 1, f, 


p f= P 1= P a 


The research reported herein was supported by 
the National Aeronautics and Space Administration 5 
under the Contract No. NAS-27799. 1 f j 

■ * K 


8 P f , p| Is the largest subset of 


that P 4 -(*:V(x) <V Hp 1 < V np b<^< Vmp ). 

4 i . 4 L 


such 


P L* P a - ? l 


l i 

{x:xeP # and . 


k t s> a b a u r 

| |P(.j boumtary of a closed set P^.j ; p , |.)* p (.)* p (.) 


; Jxl_.1-dlnn&iBDal..ruL vector. space. 
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‘-'•'X. SriHM PERIODICALS CO. 


s i " ) - [s{ -) s*‘* ... *<‘>] T » 0 , (•) • a, f, 1, ! 
prescribed constant positive vector. 

sgn x • dlag (sgn Xj sgn x 2 ... sgn x (1 ) . 

t £ R 1 , •* < tQ < t < *• , 

t: either red 1 and t > 0 , or t » , 

t s £r' , t s €(0, t) , system settling time. 

res' , t - ttstQ < t < t„ ♦ t) , t s - <tjt 0 ♦ i 

V:K n -► R 1 , V /.» • min V(x) , V H( - - max V(x) . : 

' (*) ' * (*) 
f(x) Is the total time derivative of Y(x) 
along solutions of a system. 

x • (x 1 x 2 ... x„) T , state of a system. 

x(t}Xg, t Q ): R n * T ■» R n , motion of a system sat4 
Isfylng x(t 0 ;x 0 , t Q ) ■ Xg . 

z:R n «T«lf , z(x, t) Is a disturbance vector. 

6. . Kronecker delta. 

1J ! j 

l | 

■ {+1 or :-l) - const., 1 • 1, 2, ... , n . 

t ■ dlag {tJ t 2 ... s n ) . j 

„ T ! 

o:R n « li* « |T ♦ R n , o ■ (oj » 2 ... o n ) . 

1*1 < ** 1*^1 < sj’ 1 , 1 * 1, 2, ... , n . , 

l |(*)|| • C (-) T (*)] 1/2 . h 


It Is Important to note the difference between the 
I definition of contractive stability of ! reference 

(I) and the definition proposed above. The former 
;1s related to practical stability over a finite 
time Interval and only the existence of a number 
!ij Is required. The latter Is concerned with 
practical stability over a given time interval T , 
‘which may be either finite or infinite. Moreover, 
'the nunfear t s £(0, i) Is a specified positive 
jmmfcer. In the sequel, x ■ ■*- . 

Now we state: 

Theowi 1 : Let V(x) « b^|x| . System- (1). it pane* 
; tically stable with the settling time : if these 

j exists a constant vectas b » 0 such that the iol- 
lading conditions one satisfied 

T , j 

!(i) b's 1 < min {U.s?} , i 
“J-l,2,....n | 

! 1 t 

(II) i(x) < 0 on P L « p « T , | 

(III) t $ > [V^ - V H( p^. p 1 jj v H ^ „ J „ r < 


j Theorem 1 Is proved In Appendix 1. For a class of 
I system (1) satisfying assumptions defined In the ' 
{sequel, e simple sufficient condition for the axis- 

i ! 

| tance of a vector b » 0 , which satisfies the 
| condition (11) of Theorem 1, Is presented In 
j Theorem 2. j j 


3. PRACTICAL STABILITY WITH SPECIFIED 


I 

SETTLING TIME 
of systems I 


In this pap^r, we shall study a class of systeas ! S 
governed by the vector differential equation i I 

■ b(*» *» t) « 0 j £ 

where z ■ z(x, t)eP on P # x T . Motions ; j 
x(t; Xg, t Q ) of system (1) are required to satisfy 
x(tj Xg, tg)€P # , Ytel , and x(t; Xg, tg)eP f , | 
VtcT s , whenever Xg£ P^ , zeP . More precisely:! i 

Ueilnitloni 1 System (1) is psacticaJUy stable mithl | 
the settling time t $ ii and only H XgeP ( , zeP j 

dT* x(t:xg, t 0 )£P a , ¥teT , and ) 

(ID x(t;x^, t 0 )£P f , ¥t£T s . j I 

u_ 


: A.) The vector function g(x, z, t) may be 

written In the form g(x, z, t) » H(x, z, t)*; 
f(x, z, t) , l.e., system (1) mayi be described 
by equation ’ 

JT * H(*, *• t) f(x, Z, t) •; (2) 

j A. 2 h <t (x, z, t) < 0 on P^ * V * T j, 1*l,2,...,p. 
(If this assumption Is satisfied {only for 
1*1,2,... , r < n then two additional 
j estimations should be satisfied (43). 

| A. 3 |h^(x, z, t)| « •*- on P # * P 4 T , ! 

1» J ■ 1* 2, ... , n . i 

; A.4 fy(x, z, t) sgn Xj » 0 , x^ f 0 , on 

P^ « P x T , 1 » 1, 2 n . 

! Under the assumptions' A. 1-4 a matrix A*(a^j) Is 
..defined by ■ I 
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*u ’ * 4 n 1nf l h n lx * *• *M + 
J J p l *v*t 1J 


11-4.,) sup |h.,U, z, t)| , (3); 
P. *PxT 


f. h to i> 


jji 
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138 _ 


7^° 


-69 


oj >_ 3.45 
|o, I < 3.45 


Oj _< -3.45 


, NjCoj)* 60 sgn o 2 . 


and we have: 


Theoaam 2 : li XGH thejie. ixiit inffiniAz number 

o{ vcctoAS b > 0 6uch that 

V(x) ■ b T (sgn x)x < |f T (x, z, t)|A T b < 0 on P. 
Theorem 2 Is proved In Appendix Z. 

The following procedure for an analysis of practl- 1 
cal stability with specified settling time results! 
from the previous theorems. A matrix function 
H(x, z, t) and a vector function f(x ( z, t) 
should be selected to satisfy the assumptions 
A. 1-4. Knowing H(x, z, t) and f(x, z, t) a 
matrix A * (a^j) has to be determined according 
to equation (3). The procedure may be continued j 
only If ACN , which guarantees (Theorem 2) the j 
existence of Infinite number of vectors b > 0 
such that the condition (11) of Theorem 1 Is satis} 
fled. A vector b > (1 should be chosen to satis- j 
fy the condition (1) of Theorem 1 and the Inequal-: 
Hy 


» i 

5 M°3* °3* 


10 sgn o 3 


§ i 

li 

§ ? 

J ! 

9 : 


/1« 3 I> or ; 

*ll« 3 l > 2. o 3 o 3 < 0 J 

~ „ - 4 * °3 i 3; > 

W* 2 ; J 


n 4 (o 4 ) 


° 4 I 


12 sgn , 

3(|o 4 |- 2) sgn o 4 , 2 <_ |o 4 | <.6 

0 . I» 4 I«* 



. s 1 » 

C 4 \ 


30 

5 , s f - 0.5 

70 > 


k ] °/ 

m.o ; 


, v * 12 sec. : 


f In tills example, f(x, z, t) and H(x,| z, t) are 
j Selected as f(x, z, t) » x 


A*b < 0 


(4); 


!H(x, z, t) 


j -0. 2-0.1 ^° -- ; 0 i -25o V° ~- sin 3^ 

.1 Ulilli .T-20 3 --- 


°\\ 

, I 

°z 

M 

°3 


\o 4 1 

1 \ 




Provlded that such a vector b > 0 has been deter- 
mined, the condition (111) of the same theorem has; 
to be tested. | 

4. EXAMPLE 

The outlined procedure for analysis of practical j 
stability with specified settling time Is now 
applied to 4 ys tern (1), where 

/•0.Z 0 0 

g(x, z, t) ■ ( 0 -1.0 2.0 

.0 0 -5.0/ 


/ W 

Ng (* 2 ) 
N 3 (o 3 , <? 3 ) 

W 

0 
0 
0 

10' 3 M 



j-1 h 4 (o 4 ) I 

V> i * I / 

The above selection of functions f(x, jz, t) and 
H(x, z, t) Is made In order to satisfy assumptions 
A. 1-4. According to Equation (3), we get i 

( -0.5 2 o\ ! 

0 -41 2 J C H . j 

, 5 30 - 5 / : 

;The vector b is chosen to be b ■ (1(^1 22 10) T , 

| such that 1st and 2nd condition (Theorem 2) of 
I Theorem 1 are satisfied. j 

1 i ! 

| The following nuatiers are calculated ffr the selected 
I ! vector b as i 


V M(P r P[) 


• 614 V mpte • 10 V^l ■ -50.5 . 


■Therefore, 
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tV nP* * V H(P r P|) 3 *1 


LkVxT 


< 12 


r 


Hence, the conditions (1-111) of Theorem 1 are sat-i 
isflad, and the analyzed system Is practically 
stabla with the settling time 12 sec. I 

5. CONCLUSIONS 

In this paper, a definition of practical stability 
with specified settling time Is proposed. The out- 
lined analysis provides Information about quality i 
of the forced system dynamic behavior over a pre- 1 
scribed time Interval, which may be finite or In- ; 
finite, and about a value of the system settling 
time. 

An analysis of the proposed type of practical sta- j 
blllty has been carried out for a broad class of ! 
nonlinear nonautonomous dynamic systems. The clast 
of systems Is defined by some general requirements I 
and Is limited neither by the order of system nor 
by the form of nonl Inearl ties. The latter fact 
shows that the proposed procedure enables an In- 
vestigation of "absolute” practical stability. 

The stability test Is reduced to verification of i 
the Hurwltz property of a constant matrix and to j 
a choice of a positive vector to satisfy given In- j 
equalities. In complex situations procedure mey ! 
require machine computations. 
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1 APPENDIX 1 


(A-2) 


| Therefore 'XjttjT^YTjpl hJt ttl ?j«h*T~, 

j J ■ 1, 2, ... , n , l.e. I 

* j 

| x(ts Xq, t 0 )€P # , on P^DxT . , 
ifurthermore 

| Wx)| Pl <»»<o. 

1 

; Integrating the last expression In (A-2) from t Q 
j to t , we get 

| | V(x)|p L - W(V|p r p1 i»HP t (t-y < 0.{A-3> 

I iSInce p£ Is defined as the largest of P f such 
j [that 

\ : V(x)|p1 , 


(A-4) 


[than from (A-3) It follows that 


v M(p r p{)- *| < A - s > 


‘Ixltj Xj,, tj,) ■ *0 ' T sa - 

T sa ’ (t |x(t; Xq. ty) P 4 ' t 0 ) 


where 


p r VxT s • i 


If b > 0 then the following Is true: j 

V(x) • b T |x| > 0 , Vx p 0 ; V(x) •* - for ||x||- 

and since V(x) < 0 on P^xDxT , according to con-j 
dltlon (11) of Theorem 1, then V(x) Is a system I 
Liapunov function on P^UxT . From the same con- j 
dltlon It follows that, for ¥teT , 

Vlx(ts Xy, ty)] iV(Xg) . From condition (1) of j 
Theorem 1, and from the previous result, one ob- 
tains i 

Vlx(t; Xy, t^)] iV(Xy) ^.bV * min 

J“1 ,2,...,n ! 

on P^I/xT . . . . (A-U. 


actual settling time, which Is evidently less than! 
the given t s . From (A-2) and (A-4) ft also 
follows that |Xj(t; Xy, ty)| «. sj onj e - U - T 

j B If 2f ••• g R g 1 

x(t| Xy, tp)ePf on P^VxTj 
This proves Theorem 1. 

APPENDIX 2 

The proof of Theorem 2 Is based on Perildskll's 
lemma In (3] and the following assumption: 

J : A Mumotion (A-A) : The elements of a certain basis. 

3 j (c,) and the coefficients of the matrix A • (a,.) 

I i 1 i 1 J 

I ■ are related by 


(A-6) 


Is the 


i 


j 1 j — ® • 1 ^ J , 1, J " 1 1 


2, . « . , n. 
(A-l) 

! Lama (3] : Assume the coefficient* ofj the maVUx 
i A • (a^j) ate such that the Anequatitjie* (A-l) 

I ate satisfied ion a contain basts Uj)j. Then in 
j | ojuUk ion ail toots of the i teut on equation 

det(A - xl) ■ 0 
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to have, negative neat paxu it ii neceuaty am 1 : j 

sufficient 44 (o\ any positive vectcA c a vectoA i 3 
b detemnined 4/iom the equation 


pJtb 


-cc 


! t 

(a- 2) j 

■ t 


•ci positive. 

Let the elements a^ of the matrix A be defined 
by 


*U 


■ -a 


U 


1n f _ |h^j(x ( 2 , t)| + 


P^xVkT 


; 5 

I j 

I I 


(1-4^) sup |h 1 .(x, 2 , t)| , (A-3)! 

P L *y»T 

provided that h^(x, 2 , t) satisfy the assunp- 
tlons A. 1-3, 1, J ■ 1, 2, ... , n . It can be 
easily verified that all elements a^j , Equation i | 
(A-3), satisfy Assunptlon (A-A) whenever c Is ] j 
chosen to be the Indentlty matrix, l.e. 1 ? 

s • d1ag(l 1 ... 1) * I . (A-4)j j 

Now, we proceed to prove Theorem 2. For system 
(1) satisfying the assumptions A.l-4 and for 

V(x) • b T |x| • b T (sgn x)x , 

V(x) ■ b T (sgn x)x • b T (sgn x) H(x, 2 , t) f(x,z,tM 
■ f T (x, 2 , t) H T (x, 2 , t)(sgn x)b ,(A-5)j 

Using assumption A.4 and Equations (A-3) we obtain: 
the following relation for V(x) 

• T T * 

V(x) ■ |f (x, 2 , t) | (sgn x) r(x, 2 , t)(sgn x)b < 

< |f T (x, z, t) I A T b on P u xVxT . (A-6) 

If A £N , according to the Lemma [3] and Equations 
(A-3), (A-4), for any c > 0 , there exists a vec- 
tor b » 0 such that 

A T b - -e < 0 . (A-7)j 

From Equations (A-6), (A-7) we conclude that | 

• T I 

V (x) <-|f‘(x, 2 , t)|c < 0 on P L «Pxl 

and since c > 0 Is an arbitrary positive vector^ : 
the proof of Theorem 2 Is complete. 
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